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Abstract Recent results from the KASCADE extensive air 
shower experiment are presented. After briefly reviewing the 
status of the experiment we report on tests of hadronic inter- 
action models and emphasize the progress being made in un- 
derstanding the properties and origin of the knee at Eknee — 
4 • IQis eV. 

Analysing the muon- and hadron trigger rates in the KAS- 
CADE calorimeter as well as the global properties of high 
energy hadrons in the shower core leads us to conclude that 
QGSJET still provides the best overall description of EAS 
data, being superior to DPMJET II-5 and NEXus 2, for ex- 
ample. 

Performing high statistics CORSIKA simulations and apply- 
ing sophisticated unfolding techniques to the electron and 
muon shower size distributions, we are able to successfully 
deconvolute the all-particle energy spectrum into energy spec- 
tra of 4 individual primary mass groups (p. He, C, Fe). Each 
of these preliminary energy distributions exhibits a knee like 
structure with a change of their knee positions suggesting a 
constant rigidity of i? = 2-3 PV. 



1 Introduction 

The origin and acceleration mechanism of ultra-high energy 
cosmic rays {E > 10'^'* eV) have been subject to debate for 
several decades. Mainly for reasons of the power required to 
maintain the observed cosmic ray energy density of Eci ~ 1 
eV/cm^, the dominant acceleration sites are generally be- 
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lieved to be supernova remnants (SNR). Charged particles 
mainly originating from the surrounding interstellar medium 
of the pre-supernova star may get trapped at the highly super- 
sonic shock wave generated by the SN explosion. Repeatedly 
reflections on both sides of the shock front lead to an accel- 
eration by the so-called 'first order Fermi mechanism' . Nat- 
urally, this leads to a power law spectrum dJ / dE cx E^^ 
as is observed experimentally. Simple dimensional estimates 
show that this process is limited to E'max ^ ^ x (p x B), with 
Z being the atomic number of the cosmic ray (CR) isotope 
and p, B the size and magnetic field strength of the accel- 
eration region. A more detailed examination of the astro- 
physical parameters suggests an upper limit of acceleration 
ofE^^^ w Z X IQis eV (Drury, 1994; Berezhko and Kseno- 
fontov, 1999). Curiously, the CR spectrum steepens from 
7 ~ 2.75 to ~ 3.1 at £; ~ 4 X IQi'^ eV which is called the 
'knee'. The coincidence thus may indicate that the 'knee' is 
related to the upper limit of acceleration. 

Alternative interpretations of the knee discuss a change in 
the propagation of CRs from their sources to the solar sys- 
tem. Such kind of propagation effects are conveniently de- 
scribed by an 'escape time' from our galaxy, t^sc- Extrap- 
olation of Tcsc from direct measurements in the GeV range 
to higher energies breaks down ai E k, i ■ IQ^^ eV, because 
cTesc ~ 300 pc which is the thickness of the galactic disk 
(Gaisser, 2000). The value corresponds to just one cross- 
ing of the disk and would give rise to significant anisotropics 
with respect to the galactic plane when approaching this value. 
Similarly as to the process of acceleration at SNR shocks, the 
process of galactic containment is closely related to magnetic 
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field confinement, i.e. in addition to anisotropies one again 
expects E^l^ (x Z. 

A picture related to both of these interpretations has been 
proposed by Erlykin and Wolfendale (1997). They consider 
the knee as a superposition of a weakly energy dependent 
galactic modulation with additional prominent structures in 
the flux spectrum caused by a single near-by object. This so- 
called 'single source model' assumes that a shock wave of a 
recent nearby supernova which exploded some 10,000 years 
ago at a distance of a few hundred parsecs, currently prop- 
agates (or has recently propagated) through the solar sys- 
tem causing distinct peaks of elemental groups in the en- 
ergy spectrum. The most recent update of this model can be 
found in the proceedings to this conference (pi 804). How- 
ever, there is some controversy whether or not the statistics of 
presently existing data gives sufficient support to the model. 

Besides those kind of astrophysical interpretations, also 
particle physics motivated pictures of explaining the knee 
were put forward. For example, Wigmans (2000) suggested 
the inverse /J-decay reaction p + i7e — > n -|- e+ with mas- 
sive relic neutrinos could destroy protons. Simple kinematics 
shows that this channel is open for _Ep > 1.7-10^^eV/m^(eV) 
Thus, a knee energy -Eknee — 4 PeV would correspond to an 
electron neutrino mass of my ^ 0.4 eV, a value presently not 
excluded by any other observation or experiment. However, 
'eating' sufficiently large amounts of protons by such a pro- 
cess requires extraordinary high local densities of relic neu- 
trinos, which appears doubtful even if possible gravitational 
trapping is considered. In case of a non-zero magnetic dipole 
moment of massive neutrinos, also electromagnetic interac- 
tions with much larger cross section would be possible, thus 
giving rise to inelastic GZK-Uke pz^-interactions. Such a pic- 
ture requiring neutrinos masses of about 100 eV has recently 
been discussed by Dova, Epele, and Swain (2001). 

Nikolsky (1995) suggested that the knee is not a property 
of the primary energy spectrum itself, but may be caused by 
changing high-energy interactions in the Earth's atmosphere. 
Producing a new type of a heavy particle in the first inter- 
actions escaping unseen by air shower experiments, or an 
abrupt increase in the multiplicity of produced particles (see 
proceedings to this conference, pl389) could, in principle, 
mimic a break in the spectrum. From the particle physics 
point of view this is not completely ruled out as the centre- 
of-mass (cms) energy available at the knee is above Tevatron 
energies. However, a well worked out model simultaneously 
describing accelerator and CR data without violating general 
physics, hke the unitarity principle, is still to be presented. 

Very recently, exotic scenarios with extra dimensions and 
TeV-scale quantum gravity have been discussed by Kazanas 
and Nicolaidis (2001) and others. Again, the argument is 
that the cms-energy of the knee (i/i = 2 TeV) is just above 
the highest energies reached by present accelerators and that 
some part of the CR energies is transferred into invisible 
channels, in this case gravitational energy. Some advantage 
of this picture is that the required fast growth of cross section 
for graviton radiation does not necessarily violate S-wave 
unitarity. 



A common feature of the particle physics interpretations 
of the knee is the expectation of seeing the break of different 
elements to be displaced by their mass number A rather than 
by their nuclear charge Z. This is understood from the reac- 
tion mechanism being governed, at sufficiently high energy, 
by the energy per nucleon E/A of the incident particle. 

To distinguish between these various ideas, better mea- 
surements of the energy spectrum and elemental composition 
of cosmic rays in the knee region are essential. Most com- 
monly, the composition is expressed in terms of the mean 
logarithmic mass, {In A). Obviously, such an average quan- 
tity appears to be rather insensitive to details of the individ- 
ual energy spectra. Clearly, the most sensitive testbed for 
comparing different models would be given by CR energy 
spectra of different elemental groups. Direct measurements 
could in principle provide such information. However, due 
to their very linuted collection area and exposure time, direct 
measurements across the knee energy are practically impos- 
sible. On the other hand, the well known difficulties in the 
interpretation of EAS data have not yet permitted such mea- 
surements either. Thus, beyond the knee energy very littie is 
known about CRs other than their all-particle energy spec- 
trum (Watson, 1998). 

The primary goal of the KASCADE experiment is to pro- 
vide such information by measuring the electromagnetic, mu- 
onic, and hadronic EAS components with high precision in 
each individual event and by applying advanced methods of 
multi-parameter data analyses. Determining the energy spec- 
trum and composition of cosmic rays in the knee region with 
high precision also calls for detailed EAS simulations, as 
performed by the CORSIKA package (Heck, 1998). Testing 
and improving hadronic interaction models entering the EAS 
simulations is another objective of KASCADE. As discussed 
below, this is realised mostly by investigating the properties 
of high energy hadrons in the shower core. Reconstructing 
properties of the primary CRs from different sets of EAS ob- 
servables and comparing their results provides another pow- 
erful means to judge the rehabihty of measurements and sim- 
ulations. 

2 Experimental 

KASCADE (Karlsruhe Shower Core and Array Detector) is 
located at the laboratory site of Forschungszentrum Karls- 
ruhe, Germany (at 8° E, 49° N, 110 m a.s.l.). In brief, it 
consists of three major components (see Fig. 1); 

1. A scintillator array comprising 252 detector stations of 
electron and muon counters arranged on a grid of 200 x 
200 m^ and providing in total about 500 m^ of e /j- and 
620 of /x-detector coverage. The detection thresh- 
olds for vertical incidence are Eg > 5 MeV and Efj, > 
230 MeV. 

2. A central detector system (320 m^) consisting of a highly- 
segmented hadronic calorimeter read out by 44,000 chan- 
nels of warm liquid ionization chambers distributed over 
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Fig. 1. Schematic layout of tlie KASCADE experiment (left), with its streamer tube tracking system (top right) and central detector (bottom 
right). 



9 read-out layers, another set of scintillation counters 
above the shielding (top cluster), a trigger plane of scin- 
tillation counters in the third iron gap and, at the very 
bottom, 2 layers of positional sensitive MWPC's, and a 
streamer tube layer with pad read-out for investigation 
of the muon component at > 2.4 GeV. 

3. A 48 X 5.4 m^ tunnel housing three horizontal and a two 
vertical layers of positional sensitive limited streamer 
tubes for muon tracking at iJ^ > 0.8 GeV. 

More details about the experiment can be found in Doll 
et al. (1990) and in Klages et al. (1997). First correlated 
data have been taken with some parts of the experiment since 
1996 and with its full set-up since 2000. At present, more 
than 500 Mio. events have been collected in a very stable 
mode and with a trigger threshold of the array corresponding 
toE^A- IQi^ eV. 

3 Tests of High-Energy Hadronic Interaction Models 

The observation of EAS provides an opportunity to study 
global properties of hadronic interactions in an energy range 
not accessible to man-made accelerators. For example, the 
cms-energy at the Tevatron collider corresponds to a fixed 
target energy in the nucleon-nucleon system of Ep ~ 1.7 • 
10^'^ eV. Even more importantly, the diffractive particle pro- 
duction dominating the energy flux in the forward region in- 
fluences the EAS development most strongly, but has been 
studied experimentally only at comparatively low energies of 



y/s ~ 10 GeV (Kaidalov, 1979). Most of the beam energy in 
present collider experiments remains unobserved. For exam- 
ple, the UA5 experiment could register up to 30 % of the total 
collision energy at y/s — 0.9 TeV, while the CDF detector 
registers only about 5 % at y/s = 1.8 TeV. Hence, hadronic 
interaction models applied to higher cms-energies and to par- 
ticle production in the very forward region rely on extrapola- 
tions and may cause systematic uncertainties in simulations 
of EAS. Additional uncertainties arise from simulations of 
p-nucleus and nucleus-nucleus collisions including a possi- 
ble formation of a quark-gluon plasma. Again, such data are 
important for EAS interpretations but have been studied only 
at low energies in the past (SPS and ISR at CERN). Only 
very recently, RHIC data at ^/s = 200 GeV have become 
available and will be very helpful in this respect. 

In KASCADE, hadronic interaction models have been tested 
by EAS data employing two basic approaches: At lower pri- 
mary CR energies, where data from direct measurements still 
exist, the energy spectra and rate of single hadrons not ac- 
companied by air showers were studied. At higher CR en- 
ergies, on the other hand, various features of the hadronic 
shower cores were investigated. In both approaches, the sen- 
sitivity to details of the hadronic interaction models arises 
mostly from the high energy threshold applied to the hadron 
reconstruction in the calorimeter. 

3.1 Tests at Low Primary Energies: Study of Trigger Rates 

An important source of uncertainty in hadronic interaction 
models originates from the energy dependent inelastic cross- 
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Fig. 2. Trigger rate vs hadron rate in the KASCADE central detec- 
tor. The top panel compares different hadronic interaction models 
to the experimental data. The systematic uncertainty, mostly given 
by the absolute flux uncertainty of the direct experiments, is indi- 
cated by the dotted line. The lower panel shows results obtained 
from the QGSJET model with modified inelastic cross section and 
diffraction dissociation (Risse, 2000; Antoni et ah, 2001). 



sections (Block, Halzen, Stanev, 2000). Taking into account 
the deviations among different experiments at Tevatron en- 
ergies, the proton-antiproton cross section is not known to 
better than 5 % at best (Avila et ai, 1999). As shown be- 
low, these uncertainties are amplified when predicting abso- 
lute hadron fluxes at sea-level by means of EAS simulations. 
Thus, one may take advantage of this strong dependence and 
perform stringent tests of models on the basis of EAS data. 
The idea is to use data of absolute CR fluxes up to several 
TeV of energy, as obtained from balloon- and satellite borne 
experiments at the top of the atmosphere. Taking into ac- 
count the measured energy distribution and chemical compo- 
sition, these particles are then propagated through the atmo- 
sphere employing the CORSIKA simulation package. At the 
level of the KASCADE experiment, we then ask for triggers 
released by either high energy hadrons (Eh 90 GeV) or 
a minimum number of 9 muons (E"^ ^ 0.49 GeV) detected 
in the central detector of KASCADE. Such an energy cut ap- 
plied to the hadrons automatically selects particles originat- 
ing from the first high-energy interactions in the atmosphere. 
Employing different hadronic interaction models, these sim- 
ulated inclusive trigger rates are then compared to actual ex- 
perimental data (Antoni et ai, 2001). Figure 2 shows the 
results. None of the predictions agrees well with the exper- 
imental data, particularly the hadron rates are overestimated 
by up to a factor of 3. Furthermore, there are also large dif- 
ferences found among the models. This convincingly proves 




Fig. 3. CORSIKA simulations of Feynman-x distributions of lead- 
ing baryons from p + N interactions at 10^'' eV energy. Five dif- 
ferent interaction models (versions of 1997) are shown (Heck et ai, 
2001). 



the sensitivity of the experimental observable to details of the 
interaction models. 

Earlier investigations of high-energy hadrons observed in 
the shower core (Antoni et ai, 1999) lead us to the conclu- 
sion that QGSJET (Kalmykov and Ostapchenko, 1993) pro- 
vides the best overall prescription of EAS data. This was in- 
dependently concluded also from a consistency analysis in- 
cluding data from different EAS experiments (Erlykin and 
Wolfendale, 1998). Therefore, several modifications were 
applied to the QGSJET model in order to study their influ- 
ence to the predicted trigger rates. Results are presented in 
the lower panel of Fig. 2. Increasing the inelastic proton-air 
cross section by 5 % (10 %) reduces the predicted hadron rate 
by approx. 27 % (54 %). Similar effects are obtained by low- 
ering the diffraction dissociation by up to 10 % of the inelas- 
tic cross section. Clearly, there is a strong sensitivity to these 
parameters. Since the total inelastic cross section appears 
to be the better known quantity, we conclude that the uncer- 
tainty arises mostly from the diffraction dissociation which 
may be overestimated in the simulated hadron-nucleus inter- 
actions by about 5 % of the inelastic cross section crf^^J^" ■ In 
fact, the trigger rates predicted by the different models shown 
in Fig. 2 (top) are strongly correlated to the simulated shapes 
of the Feynman-a; distributions of leading baryons. This is 
demonstrated in Fig. 3 for CORSIKA simulations of p + 
interactions at £' = 10^® eV; for the same total cross sec- 
tion, a large fraction of diffractive interactions leads to high 
trigger rates and vice versa. In order to verify the aforemen- 
tioned conclusions about the observed deviations of the trig- 
ger rates it would be highly desirable to measure the proton 
inelastic cross section on nitrogen and oxygen at the highest 
energies at accelerators and to study the kinematical region 
of diffraction dissociation. 

3.2 Tests at High Primary Energies: Study of Hadrons in 
the Shower Core 

While the previous study selects mostly diffractive type in- 
teractions in the atmosphere, the following analysis samples 
the more 'typical' nuclear interactions leading to a fully de- 
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Fig. 4. Deviations of the simulated hadronic energy sum normal- 
ized to experimental KASCADE data for QGSJET, NEXus 2, and 
DPMJET 11-5. Data in the upper and lower panel are binned as 
function of electron- and truncated muon size, respectively. (Milke, 
2002) 



tation best. With some exceptions at lower electron sizes, 
also DPMJET provides a reasonable description of the ex- 
perimental data. However, NEXus 2 clearly fails in this test: 
there is too little hadronic energy at given electron size. Since 
hadrons in an EAS continuously populate the electromag- 
netic component, this mismatch points to a general problem 
of the balance between electromagnetic and hadronic energy 
in the used version of NEXus. However, it should be men- 
tioned that this version of NEXus is still considered prelim- 
inary and will be subject of further improvements. The ob- 
served overall decrease of simulated/experimental hadronic 
energies may be due to an increasingly heavier composition 
in the experimental data at higher energies (see next section). 

The truncated muon number provides an energy estimator 
which is almost independent from primary mass. Thus, re- 
peating the analysis and plotting the same quantities now as 
a function of truncated muon size should yield results with 
the proton and iron simulations being reasonably above and 
below the experimental data. Again, QGSJET exhibits the 
expected behaviour (lower panel of Fig. 4). Different from 
above, now the NEXus model follows the QGSJET results 
very closely and DPMJET fails the test; DPMJET II-5 pro- 
duces significantly too many hadrons so that a composition 
based on this model would appear too heavy. Such a compo- 
sition contradicts measurements from other experiments as 
well as those from KASCADE (see below). At Ig Nf^' ^ 4.5, 
i.e. energies corresponding to E ^ 10^^ eV, a composition 
even heavier than iron would be needed to describe the exper- 
imental data. The combined results of Fig. 4 top and bottom 
thus suggest that EAS simulated with DPMJET II-5 penetrate 
too deeply into the atmosphere. The conclusion is confirmed 
also by studying the multiplicity and energy distribution of 
hadrons in the calorimeter (Milke et al, 2001). 



ve loped EAS. Events have to fulfil standard cuts in electron 
and muon number for EAS (Ne > 10^, NfJ' > 10^) and need 
to have at least one hadron with E ^ 50 GeV reconstructed 
in the calorimeter Details of this type of analysis were pre- 
sented by Milke et al. (2001) at this conference. An exam- 
ple employing only the three models which are closest to the 
data in Fig. 2 is depicted in Fig. 4. Here, the total hadronic 
energy observed in the calorimeter is compared with COR- 
SIKA simulations for proton and iron primaries employing 
the QGSJET (Kalmykov and Ostapchenko, 1993), NEXus 2 
(Drescher et al., 2000), and DPMJET (version II-5) (Ranft, 
1999) models. The upper panel shows the ratio (Ig ^ E^si^- 
hY. Eh,data) / hY. EhAata for cach of thcsc modcls as a 
function of electron size. The lower panel shows the same 
quantity as function of the truncated muon size, defined as 
iV^'' = J^^^™ 2Trpfj^{r)dr. Considering the steep primary 
energy distribution and the different EAS developments for 
light and heavy primaries, binning the data as a function 
of Ig Ne tends to enrich showers originating from light pri- 
maries. Therefore, proton simulations (open symbols in the 
upper panel of Figure 4) are expected to closely resemble the 
experimental data. QGSJET simulations follow this expec- 



4 Determination of the primary energy spectra and mass 
composition 

Results presented in the previous section have demonstrated 
that QGSJET is still the model of choice when analysing 
EAS data in the knee region. Deviations between the Monte 
Carlo predictions and experimentally observed distributions 
of high energy hadrons {Eh ^ 90 GeV) have become rea- 
sonably small, however, efforts are still needed for further 
optimization of interaction models. As pointed out before, 
selecting high energy hadrons amplifies the sensitivity to the 
actual modelling of high energy hadronic interactions. Low 
energy electrons and muons, on the contrary, should be less 
affected as they result from an average of many interactions 
at mostly lower cms-energies. Thus, in the following we 
will concentrate on the reconstructed electron and truncated 
muon numbers. 

4. 1 Unfolding Techniques 

It is well known that for given zenith angle, i.e. fixed atmo- 
spheric thickness, and Nj^'^ depend simultaneously on the 



6 



KASCADE preliminary 



••••• 



■ sum of all 

• proton 

A helium 

* carbon 
T iron 




Primary Energy (eV) 



Fig. 5. Preliminary energy distributions 
of four primary mass groups as obtained 
from tiie unfolding procedure (Ulrich et 
ah, 2001). The sum of the 4 individual 
distributions represents the all-particle 
CR-spectrum and is shown by the black 
squares. The vertical error bars repre- 
sent the statistical uncertainties which 
are dominated by the statistical uncer- 
tainties of Monte Carlo simulations (see 
also Fig. 6). 




Primary Energy (eV) 



Fig. 6. Results of repeated unfolding procedures using different 
input parameters from the Monte Carlo simulations. The scattering 
of the symbols has been used to estimate the statistical uncertainty 
of the Monte Carlo simulation. See also Fig. 5 and text for details 
(Ulrich, 2002). 

mass. A, and energy, E, of the primary particle. Besides their 
specific average relation {N^, N^) <-> {E, A) also their^Mc- 
tuations are known to change with energy and mass. Clearly, 
such effects are to be taken into account when reconstruct- 
ing the properties of the primary particles, particularly in the 
presence of steeply falling energy spectra. Mathematically, 
the observed electron and muon shower size distributions at 
a given zenith angle bin can then be written as 

where the sum runs over all primary masses A. The quan- 
tity PA (Ig Ne,f_i I Ig E) denotes the probability for a primary 
particle of mass A and energy Ig E to be reconstructed as an 
air shower with electron- and muon size IgiVg and \gNjf, 



respectively. The equation thus represents a set of Fredholm 
integral equations of 1*^' kind. In practice, not all primary 
masses A will be taken into account, but only a representa- 
tive group of different primary particles, here protons, he- 
hum, carbon, and iron. Suitable methods to solve the in- 
verse equation, i.e. to infer the physical quantities of interest, 
(IJa /d\gE, from the experimentally observed dJ/d Ig A^e./i 
distributions, are unfolding algorithms. Here, we have cho- 
sen the iterative scheme of the Gold-algorithm (Gold, 1964) 
which tries to minimize the functional and allows only 
non-negative solutions (Ulrich et al., 2001). 

The kernel function pA needs to be determined from EAS 
and detector simulations and has to account for all kinds 
of physical and experimental effects, such as fluctuations of 
shower sizes, trigger and detection efficiencies or reconstruc- 
tion accuracies. It has been determined from high statistics 
CORSIKA/QGSJET simulations generated with thinning op- 
tion for 3 different bins of zenith angle (0°-18°, 18°-25.8°, 
25.9°-32.3°), four primary masses and for a large number 
of fixed primary energies. The individual and N^^^ dis- 
tributions obtained for each primary energy and mass were 
then parametrized and the variation of the parameters with 
shower size, zenith angle, primary energy, and primary mass 
was used in the analysis. Details of this method are described 
in Uhich et al. (2001) and by Ulrich (2002). 

Results based on the experimental Ig and Ig Njf shower 
distributions of all 3 zenith angle bins are presented in Fig. 5 
for 4 primary mass groups. Each of the reconstructed energy 
distributions shows a knee like structure which is shifted to- 
wards higher energies with increasing mass. The knee in the 
all-particle spectrum at about 4 PeV appears to me mostly 
caused by the break of the light elements proton and helium. 
This is well in agreement to earlier analyses of KASCADE 
data (Glasstetter et al, 1999; Antoni et al., 2002). 

It should be stressed that no assumption about spectral 
shapes or mass abundances is made in the analysis presented 
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here, only the measured electron and muon shower size dis- 
tributions at 3 zenith angles and the simulated response func- 
tion PA are used. The error bars in Fig. 5 represent statistical 
errors only and are dominated by the limited statistics of the 
presently available CORSIKA simulations. Statistical errors 
of the experimental data correspond to about 20 months of 
data taking and are very small compared to that of the simu- 
lations. 

In order to investigate the robustness of the reconstructed 
energy spectra, the unfolding algorithm has been repeated 
many times employing different sets of input parameters to 
the kernel function p^. These parameters describing the shapes 
of the simulated and iV^'' distributions for fixed primary 
energies (see above) have been randomly distributed accord- 
ing to their individual fit uncertainties assuming Gaussian er- 
ror shapes. The result of this study is shown in Fig. 6 for the 
reconstructed proton and helium energy spectra. The shoal 
of symbols, all obtained using the same experimental data, 
clearly demonstrates the presence of a knee at distinct ener- 
gies in both of the reconstructed spectra. A closer inspection 
of the spectra below the knee energy exhibits a steeper pro- 
ton spectrum as compared to helium and a smoother turnover 
into a power law behaviour above the knee. However, it 
should be kept in mind that the reconstructed proton spec- 
trum is significantly affected by the large shower fluctuations 
limiting the experimental energy resolution. In fact, simula- 
tions have shown that even a sharp break would result in a 
smoothness similar to the observed one (Ulrich, 2002). The 
observed scattering of the reconstructed fluxes at a given en- 
ergy has been used to estimate the statistical uncertainties of 
the unfolding procedure. Adding this in quadrature to the 
experimental statistical errors yields the values shown in fig- 
ure 5. Clearly, improving the statistics of the Monte Carlo 
simulations in the near future will reduce this error. 

The still preliminary energy spectra support the picture of 
a rigidity dependent scaling of the knee position. However, 
since Z/A ~ 0.5 for all shown mass groups with A > A, 
the crucial test for distinguishing between a E/A and E/Z 



scaling of the knee is given by the comparison of the proton 
and helium spectra. As can be seen from Fig. 6, the knee 
between p and He is shifted by factor of w 2.5 favouring the 
rigidity picture. However, because of the different spectral 
shapes, such a comparison is not straightforward and further 
investigations are needed to verify that conclusion. Also, the 
sensitivity to the used hadronic interaction model needs to be 
studied. This may shift the overall energy scale and flux of 
the individual curves but is anticipated to have only minor 
influence to the relative positions of the observed knees. 

4.2 Chemical Composition and Comparison to Neural Net- 
work Approaches 

Individual energy spectra of the type shown above are ide- 
ally suited for comparisons with cosmic ray source- and ac- 
celeration models. However, up to now most EAS experi- 
ments have only been able to extract the all-particle energy 
spectrum and the mean logarithmic mass, (In A). As a con- 
sequence, most theoretical papers present their results just 
in terms of these rather insensitive and inclusive variables. 
Thus, for completeness we plot in Fig. 7 the quantity {in A) 
as a function of the primary energy. As expected from the 
structure of the individual energy spectra, we find an increas- 
ingly heavier composition at energies above the knee. The 
data are compared to theoretical calculations by Erlykin and 
Wolfendale (1997), Berezhko and Ksenofontov (1999), and 
Swordy (1995). Also, other experimental results from KAS- 
CADE are included. Common to all models and data is an 
increasingly heavier mass at energies above the knee. Only 
the single source model of Erlykin and Wolfendale (1997) 
shows some structure at ^ 10^^ eV due to the assumed oc- 
currence of an iron peak. As can be seen from that figure, 
the models differ significantly in their prediction of the mean 
mass and, except for the Swordy calculation, predict heavier 
masses than extracted from the experimental data. 

In addition to the results from the unfolding procedure, we 
include earlier results from an analysis of the IgiV^/ IgiVe 
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ratio (Weber et ah, 1999) and from a neural network algo- 
rithm applied to the reconstructed electron and muon num- 
bers at 3 different zenith angels (Roth et al, 2001). Within 
their errors, we find similar results from the unfolding pro- 
cedure and from the Ig N^/ Ig Ng ratios. However, the neu- 
ral network approach yields a similar trend only at energies 
above the knee. In contrast to the other results, the compo- 
sition shows a weak decrease in the mean logarithmic mass 
by about 0.5 units in the energy range from about 1 PeV to 

4 PeV. Similar trends have been observed also by Cherenkov 
air shower measurements. (For a compilation of several data 
sets see Swordy et al. (2002)). At present, it cannot be said 
whether the effect is real or an artefact of the data analysis, 
such as due to improperly accounted EAS fluctuations, neu- 
ral network training effects, etc. . For example, it is easy to 
realize that the occurrence of a break in the individual energy 
spectra combined with the weU known stronger shower fluc- 
tuations for Ught primaries could qualitatively account for 
such an effect. Generally, this would produce a bias towards 
a Ughter composition the steeper the power law distribution 
of the true proton spectrum is. 

The experimental data match well to direct measurements 
recently reported by the JACEE and RUNJOB experiments 
up to ~ 10^^ eV (Asakimori et al, 1998; Apanasenko et 
al, 2001). While RUNJOB reports an almost constant com- 
position of (In A) c:^ 1.6 in the energy range 10^^ - 10^^ eV, 
that of JACEE tends to increase from 1.5 to 2.5 in the same 
energy range. This may point to some methodical problems 
in at least one of the two experiments. However, the statisti- 
cal and systematic uncertainties of both experiments are still 
sufficiently large at lO^^ eV ((In A) ~ LSt^ ° for RUNJOB 
and 2 . 5 ± . 6 f or JACEE) to be compatible to one another and 
to the reported KASCADE data. Even fliough only 40 % of 
the available RUNJOB data have been included in that anal- 
ysis and some more data being available from JACEE, the 
unsatisfactory situation of lacking statistics in direct experi- 
ments around 10^^ eV will hardly change in the near future. 

5 Discussion and Outlook 

We have presented results on tests of the latest versions of 
currently available hadronic interaction models implemented 
into CORSIKA. Only three of the models (QGSJET, DP- 
MJET II-5, and neXus 2) could describe the experimen- 
tally observed trigger rates at a reasonable level. The mod- 
els SibyU 1.6, VENUS, and HDPM overestimate die hadron 
trigger rates by up to a factor of 3. The differences in the 
predicted trigger rates are strongly correlated to the amount 
of diffractive dissociation, as can be seen from the Feynman- 
X distributions of leading baryons in Fig. 3. Clearly, there 
is a need to improve the knowledge about diffractive inter- 
actions up to the highest energies by conducting either ded- 
icated experiments or by increasing the forward acceptance 
of currently operating and planned accelerator experiments. 

Investigating the properties of QGSJET, neXus, and DPM- 
JET in EAS at energies around the knee reveals some prob- 



lems also for the latter two models. neXus 2 predicts too 
Uttle hadronic energy at given electron size, i.e. the balance 
between electromagnetic and hadronic energy disagrees with 
experimental data. DPMJET II-5, on the other hand, pro- 
duces too many hadrons (and electrons) at given muon size, 
i.e. the EAS penetrate too deeply into the atmosphere. Thus, 
applying DPMJET II-5 for composition studies would mimic 
a too heavy composition. Combining aU results, QGSJET 
stiU provides the best overall description of EAS in the knee 
region. Nevertheless, a continuation of such tests up to higher 
primary energies with optimized experimental observables 
and with tuned cuts to the hadron energy, appears mandatory 
in absence of adequate accelerator experiments. 

Analysing the energy spectrum and composition of CRs in 
the knee region calls for suitable analysis techniques which 
properly account for mass and energy dependent EAS fluc- 
tuations and which enable the simultaneous usage of many 
(correlated) EAS parameters. Here, we have presented re- 
sults from an 'inclusive' unfolding technique applied to the 
electron and muon shower size distributions at different zenith 
angles as well as results from an 'exclusive' neural network 
algorithm. The notations 'inclusive' and 'exclusive' refer to 
the fact that the unfolding algorithm does not aim at provid- 
ing information about the primary energy and mass in sin- 
gle events, but only for the full sample of events. In other 
words, the technique is applied to measured shower size dis- 
tributions and it supplies energy distributions. This lack of 
information in single events appears to be well compensated 
for by the sensitivity to details of the primary energy spectra. 
Employing this approach provides, for the first time, primary 
energy spectra for different mass species. Without making 
any 'a priory' assumption about the individual energy spec- 
tra, we find power-law distributions with a knee-Uke struc- 
ture in each of the distributions. The observed shift of the 
knee energies supports the picture of a rigidity dependent 
cut-off, as expected from most astrophysical models of CR 
origin and acceleration. 

The neural network approach confirms these overall trends 
but quantitative differences to the unfolding technique are 
observed and deserve further studies. Using only the electron 
and muon shower size as input does not allow to adequately 
distinguish more than two or at most three mass groups on 
a shower-by-shower basis. However, another advantage be- 
sides the shower-by- shower information is the simple inclu- 
sion of more EAS observables and the automatic consid- 
eration of correlations among the parameters. Such multi- 
parameter analyses have revealed systematic effects in the 
reconstructed energy and mass depending of the EAS ob- 
servables used (Antoni et al, 2002), thereby demonstrating 
the need for further improvements of the EAS simulations. 

Clearly, much more work is still needed to verify and man- 
ifest the important finding of the observed break in the indi- 
vidual energy spectra. Despite the enormous progress made 
in recent years, parallel efforts are still required to (i) im- 
prove the modelling of EAS and to understand also the tails 
in the distributions of EAS observables, (ii) to continue de- 
veloping advanced and complementary EAS analysis tech- 
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Fig. 8. Distribution of the muon production iieigiit, expressed in 
terms of the radial angle for 3.75 < IgA'^^'^ < 4.0. The data are 
compared to CORSIKA/QGSJET simulations of proton and iron 
primaries (Biittner et ai, 2001). 



niques, (ill) to provide additional EAS observables by the 
experimental data, and (iv) to improve the statistics and re- 
construction accuracies particularly at the highest energies 
(~ 10^^ eV) in order to verify the occurrence of an iron 
knee. All of these aspects are presently addressed by the 
KASCADE collaboration. For example, besides improve- 
ments to CORSIKA and its interaction models, a new two- 
dimensional non-parametric unfolding algorithm has been 
developed to account also for the shower-by-shower corre- 
lations of the electron and muon numbers (Roth, 2002; An- 
toni et at., 2002). The preliminary results obtained with that 
method agree well with those presented in this article. Also, 
new observables just become available in KASCADE and 
will be included in future analyses. As an example we show 
in Fig. 8 preliminary results of the muon production height 
as extracted from the data of the streamer tube tracking de- 
tector in the tunnel north of the central detector (Biittner et 
ai, 2001). This observable will become very powerful par- 
ticularly at energies ^ 10^^ eV, but further optimization of 
the tracking algorithm and cut parameters is still needed be- 
fore including the observable in the standard set of energy 
and composition parameters. Finally, in order to improve 
the measurements of primary energy and mass up to energies 
of about 10^^ eV, the KASCADE and EAS-TOP experiment 
have been combined at the site of Forschungszentrum Karls- 
ruhe to form a new experiment, called KASCADE-Grande 
(Bertaina et ai, 2001). After finishing the deployment of 37 
stations of 10 m^ detection area each over an area of 0.5 
km^ by the end 2001 the new experiment will start to take 
data in spring 2002. 
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